Ambipolar spin diffusion and D'yakonov-Perel' spin relaxation in GaAs quantum wells 
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We report theoretical and experimental studies of ambipolar spin diffusion in a semiconductor. A 
circularly polarized laser pulse is used to excite spin-polarized carriers in a GaAs multiple quantum 
well sample at 80 K. Diffusion of electron and spin densities is simultaneously measured using a 
spatially and temporally resolved pump-probe technique. Two regimes of diffusion for spin-polarized 
electrons are observed. Initially, the rate of spin diffusion is similar to that of density diffusion 
and is controlled by the ambipolar diffusion coefficient. At later times, the spin diffusion slows 
down considerably relative to the density diffusion and appears to be controlled by a non-constant 
(decreasing) spin diffusion coefficient. We suggest that the long-time behavior of the spin density 
can be understood in terms of an inhomogeneous spin relaxation rate, which grows with decreasing 
density. The behavior of the spin relaxation rate is consistent with a model of D'yakonov-Perel' 
relaxation limited by the Coulomb scattering between carriers. 

PACS numbers: 72.25.Dc, 78.47.jc, 72.25.Fe, 72.25. Rb, 73.63.Hs 



I. INTRODUCTION 



Recently, concerted efforts have been directed at de- 
veloping new ways of understanding the spin-degree of 
freedom in semiconductors and its potential use in elec- 
tronic devicesji^ Although electric spin injection has 
been achieved in a number of semiconductors, such as Si, 4 
GaAsj££ InAsjZiSi carbon nanotubesj 9 -^ and grapheneji 1 - 
spin transport in semiconductors is more conveniently 
studied by optical techniques. Optical spin injection 
and detection are very efficient due to well-established 
spin selection rules^ The transport of spin can be 
studied directly with high temporal and spatial res- 
olution. For example, spin transport in GaAs has 
been studied by measuring polarization properties of 
photoluminescencei 13 ' 14 ! 15 ! 16 ! 17 Faraday and Kerr rota- 
tion techniques have been used to study transport of 
spins injected by optical orientation ) 18 i 19 i 20 by the spin 
Hall effect ; 21 ' 22 ! 23 and from ferromagnets^ Very re- 
cently, the spin-grating technique has been used to study 
pure spin diffusion i 25 i 26 i 27 i 28 ' 29 Spin transport in semi- 
conductors has also been an active topic of theoretical 
studies in recent yeaTS 30M l &&M&£2£L28£M° 

In optical excitations, the injection of spin-polarized 
electrons is inevitably accompanied by the injection of 
an equal number of positively charged and spin-polarized 
holes^ In previous experimental studies, however, the 
effect of the holes was minimized or eliminated by various 
strategies. In studies of n-type doped semiconductors, for 
instance, the optically injected holes quickly lose their po- 
larization and recombine with majority carriers, leaving 
behind a purely electronic spin packet i 18 ' 19 i 20 i 21 i 22 i 23 i 24 
In intrinsic samples with an externally applied elec- 
tric field, holes are spatially separated from the elec- 
trons, and the transport is dominated by drift of spin- 
polarized electrons in the fiekh 13 ' 14 In spin-grating ex- 
periments, interference of two oppositely polarized laser 



beams results in polarization modulation - without in- 
tensity modulation, across the laser spot - which is typ- 
ically much bigger than the grating period. Since the 
excited carrier density is uniform, no carrier diffusion 
occurs. Spin diffusion in this configuration is not in- 
fluenced by holes i 25 ' 26 ' 27 ' 28 ' 29 Finally, in spin-transport 
experiments using surface acoustic waves, spatially mod- 
ulated piezoelectric fields of surface acoustic waves sep- 
arate and move spin-polarized carriers. Since electrons 
and holes are trapped at the minima and maxima of 
the acoustic wave, respectively, holes are separated from 
electrons over macroscopic distances ! 15 ' 16 ! 17 In all these 
cases, since the influence of the holes is suppressed, spin 
transport is dominated by the drift and diffusion of spin- 
polarized electrons, analog to unipolar carrier transport. 

In this paper, we report theoretical and experimen- 
tal studies of diffusion of optically injected spin-polarized 
carriers in an intrinsic semiconductor without an exter- 
nally applied electric field. In this case, the spin diffu- 
sion is strongly influenced by the presence of the holes. 
In the experiments, spin-polarized electron-hole pairs are 
excited in GaAs quantum wells (QWs) by a tightly fo- 
cused and circularly polarized laser pulse. Since hole 
spin relaxation is shorter than a few picoseconds) 41 - after 
a short transient process, the carrier system has three 
components, spin-up electrons, spin-down electrons, and 
holes, as illustrated in Fig. la. All three species diffuse 
in the quantum-well plane and interact with each other. 
Figure lb illustrates the carrier system at a later time. 
Although holes are unpolarized, they influence the diffu- 
sion of spin-polarized electrons via the Coulomb attrac- 
tion. Basically, the rate of diffusion is controlled by the 
slow holes, while the highly mobile electrons quickly ad- 
just to neutralize the holes. This forces both spin-up 
and spin-down electrons to diffuse in the same direction, 
in marked opposition to unipolar spin packets in which 
electrons of opposite spin orientations diffuse in opposite 
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sion, where the spin diffusion coefficient decreases signif- 
icantly with time, remains for the time being beyond the 
reach of our drift-diffusion theory and will be the subject 
of further investigation. 

This paper is organized as follows. In Sec. II, we dis- 
cuss experimental techniques used for optical injection 
and detection of spin-polarized carriers. Experimental 
results are presented in Sec. III. A drift-diffusion theory 
is outlined in Sec. IV, followed by theoretical and ex- 
perimental results on an inhomogeneous spin relaxation 
time in Sec. V. We conclude with a summary in Sec. VI. 
More details on spin diffusion matrix are discussed in the 
Appendix. 



II. EXPERIMENTAL TECHNIQUES 



FIG. 1: (Color online) Experimental geometry (a) Spin- 
polarized electrons (e) and holes (h) are injected in GaAs 
QWs by a circularly polarized pump pulse. Right after the 
hole spin relaxation, the carrier system is composed of spin-up 
electrons, spin-down electrons, and holes, within the excita- 
tion spot. The arrows on the electrons indicate the direction 
of the electron spin. This is the initial state of the ambipo- 
lar carrier and spin diffusion process to be studied, (b) The 
carrier system at a later time. Diffusion of electron-hole pairs 
is driven by the density gradient, and is temporally and spa- 
tially resolved by measuring the differential transmission of a 
linearly polarized probe pulse (c) Experimental set-up. The 
pump and probe pulses are counter propagating and are fo- 
cused to the sample by microscope objectives with numerical 
aperture of 0.4. A portion of the transmitted probe pulse is 
reflected by a beam splitter to a silicon photodiode to detect 
electron density. The spin density is detected by analyzing 
the polarization state of the other portion of the transmitted 
probe pulse by an optical bridge composed of a quarter-wave 
plate (A/4), a Wollaston prism, and a balanced detector. 



directions. 

We study this triple-polar diffusion process by a pump- 
probe technique. Dynamics of electron density and spin 
density are spatially and temporally resolved by mea- 
suring carrier-induced changes in the transmission of a 
linearly polarized probe pulse [Fig. lb]. We observe a 
sub-linear-expansion process of the area of the spin den- 
sity packet, while the simultaneously measured electron 
density packet expands linearly. This indicates that the 
spin transport cannot be described as a classical diffusion 
process with a constant diffusion coefficient. The spin dif- 
fusion is significantly slower than the ambipolar carrier 
diffusion i.e., diffusion of holes screened by electrons. Our 
theoretical analysis, based on a three-component drift- 
diffusion equation, shows that the long-time behavior of 
the spin density can be understood in terms of a spin 
relaxation rate that grows with decreasing density. This 
behavior is consistent with a model of D'yakonov-Perel' 
relaxation limited by the Coulomb scattering between the 
carriers. However, the short-time behavior of the diffu- 



The GaAs multiple-quantum-well sample contains 40 
periods of 10-nm GaAs layers sandwiched by 10-nm 
Alo.7Gao.3As barriers. During the measurements, the 
sample is cooled to 80 K. Spin-polarized carriers are ex- 
cited by a circularly polarized 250-fs pump pulse with a 
central wavelength of 1550 nm obtained from an optical 
parametric oscillator pumped by a Ti:sapphire laser at 
80 MHz. The circular polarization is achieved by a broad- 
band quarter-wave plate and a linear polarizer, with a pu- 
rity better than 97% on the sample. The pump pulse is 
tightly focused to excite carriers with a Gaussian spatial 
profile of wq—1.5 /im (full width at half maxima) through 
two-photon absorption, with an excitation excess energy 
of about 40 meV. The density of spin-up electrons ex- 
cited, n-f, is approximately three times higher than the 
density of spin-down electrons, n^, owing to selection 
rules for two-photon absorption of circularly polarized 
lights The advantage of using two-photon absorption 
process instead of one-photon absorption to excite car- 
riers is that, due to a small two-photon absorption co- 
efficient, the peak electron density is uniform across the 
quantum wells since the intensity of the pump pulse does 
not attenuate significantly as it propagates through the 
sample. 

The carrier and spin dynamics are monitored by simul- 
taneously measuring electron and spin density profiles as 
functions of time by using a spatially and temporally re- 
solved pump-probe technique. The total electron density, 
n = n-j + rii, is measured by focusing a linearly polar- 
ized 100-fs probe pulse on the sample with a spot size 
of 1.5 /im [Fig. lc]. The probe pulse is obtained from 
the Thsapphire laser that is used to pump the optical 
parametric oscillator, and is tuned to Is heavy-hole ex- 
citonic resonance (803 nm) to efficiently probe carriers 
based on well-established excitonic absorption saturation 
caused by free carriers.— The differential transmission, 
AT/To ee [T{n) - T ]/T , i.e., the normalized difference 
between transmission in the presence of carriers [T(n)] 
and without them [To], is measured by reflecting a por- 
tion of the transmitted probe pulse to a silicon photodi- 
ode [see Fig. lc]. This quantity is directly proportional 
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to the electron density under our experimental condi- 
tions (see discussion below). A lock-in amplifier is used 
to detect the signal with the pump intensity modulated 
at 821 Hz by a mechanical chopper. 

Spin density, s = n-f — nj,, is measured by analyz- 
ing carrier-induced circular dichroism of the same probe 
pulse, i.e., the absorption difference of right- and left- 
hand circularly polarized probe pulses in the presence 
of spin-polarized carriers. The linearly-polarized probe 
pulse is composed of two circular components. Due to 
spin-selection rules, each component preferentially senses 
one spin system^ By using a quarter- wave plate (A/4), 
the two circular components are converted to two or- 
thogonal linear polarizations, which arc then separated 
by a Wollaston prism and sent to two photodiodes of 
a balanced detector. The output of the balanced de- 
tector is proportional to the difference between the dif- 
ferential transmissions of the two circular components, 
(AT + — AT~)/Tq, which is proportional to the spin den- 
sity 544^5,46^7 

In order to deduce the actual electron and spin den- 
sities from the measured differential transmission signals 
AT /T and (AT+ - AT~)/T , separate measurements 
are performed with the same set-up, but with a pump 
pulse of a central wavelength at 775 nm, which is ob- 
tained by the second-harmonic generation of the 1550- nm 
pulse. This pump pulse excites carriers via one-photon 
absorption. The linear relation between AT /To and n is 
verified by measuring AT /To as a function of the energy 
fluence of this pump pulse, which is proportional to the 
electron density in the linear regime. Furthermore, we 
measure the absorption coefficient of the sample at this 
pump wavelength by comparing the incident, reflected, 
and transmitted powers of the pump pulse. This allows 
us to determine the number of photons absorbed by the 
sample for a given pump energy fluence. Since other ab- 
sorption processes are much weaker than the interband 
absorption, approximately the same number of carriers 
is excited. Therefore, we obtain the proportionality con- 
stant between AT '/To and n. The constant of proportion- 
ality between the spin density s and (AT + — AT~)/To is 
similarly deduced, using in addition the well-established 
fact that the one-photon absorption of circularly polar- 
ized light tuned close to the band gap of GaAs excites 
electrons with a spin polarization of s/n = 0.5 J£ 



III. EXPERIMENTAL RESULTS 

With the experimental techniques discussed in Sec. II, 
we simultaneously measure electron density and spin den- 
sity as functions of space by scanning the probe spot 
across the pump spot, and as functions of time by scan- 
ning the time-delay between the probe and the pump 
pulses. Figure 2 summarizes the dynamics of electron 
density and spin density measured with a peak electron 
density of 2.3 x 10 17 cm -3 . The contour map in Fig. 2a 
shows how electron density varies with time and space 
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FIG. 2: (Color online) Spatiotemporal dynamics of electron 
(n) and spin (n) densities measured at a sample tempera- 
ture of 80 K and a peak electron density of 2.3 x 10 17 cm -3 , 
(a) Electron density as functions of time and space. The 
electron density is deduced from the measured AT/Tq. (b) 
Spatial profiles of electron density measured at 5 ps (squares) 
and 300 ps (circles). The profiles are normalized in order to 
show the expansion caused by the carrier diffusion. The solid 
lines are Gaussian fits, (c) Squared width of electron den- 
sity profiles as a function of time obtained by Gaussian fit 
to profiles measured at various times. The linear fit (solid 
line) corresponds to an ambipolar carrier diffusion coefficient 
of D a = 21 cm 2 s _1 . (d)-(f) Spatiotemporal dynamics of spin 
density obtained in the same scans. 



after injection. Here r is defined as the distance between 
the centers of the pump and probe spots, and t is the 
time delay of the probe pulse with respect to the pump 
pulse. Two normalized spatial profiles of electron den- 
sity measured at t=5 ps (squares) and 300 ps (circles) 
are plotted in Fig. 2b. The solid lines are Gaussian fits 
to the data. Clearly, the spatial profile of electron density 
remains Gaussian shape, and expands due to the trans- 
port. To quantitatively describe the transport, we deduce 
the time variation in the full width at half maxima of the 
electron density profile, w n , by fitting the electron den- 
sity profiles measured at all probe delays. In Fig 2c, the 
squared width as a function of t is plotted. 

According to the classical diffusion model with a con- 
stant diffusion coefficient, the squared width increases 
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linearly as 



<{t) 



161n(2)D„i, 



(1) 



where D a is the ambipolar carrier diffusion 
coef5cient i 48 i 49 By linear fit to the data, we de- 
duce D a = 21 cm 2 s . It is worth mentioning that, in 
the experiment, due to a finite size of the probe spot, 
the measured profiles shown in Fig. 2b are actually 
convolutions of the probe spot and the actual electron 
density profiles. However, since both the probe spot 
and the electron density profiles are Gaussian, the 
convolution does not influence the measurement of D a £Q 

The right half of Fig. 2 shows the spin component of 
the diffusion process measured simultaneously with den- 
sity. The contour map in Fig. 2d shows the spatiotem- 
poral dynamics of spin density after optical excitation. 
Spin diffusion is evident by comparing the two profiles 
in Fig. 2e measured at t=5 ps (squares) and 300 ps (cir- 
cles). The profile remains Gaussian, as confirmed by the 
fits (solid lines). However, the spin density profile at 
i=300 ps is narrower than the electron density profile 
shown in Fig. 2b. This indicates the different transport 
behavior of spin and electron densities. Quantitatively, 
the squared width of the spin density profile does not in- 
crease linearly, as shown in Fig. 2f, in striking contrast to 
the expansion of the electron density profile obtained in 
the same measurement. Initially, the expansion rate of 
spin density profile is similar to that of electron density 
profile. At later times, the spin diffusion slows down con- 
siderably relative to the density diffusion. The observed 
sub-linear expansion of spin density profile shows that 
spin diffusion cannot be described as a classical diffusion 
process with a constant diffusion coefficient. 

The procedure summarized in Fig. 2 is used to sys- 
tematically study the influence of electron density on the 
diffusion process. The sub-linear expansion of the spin 
density profile is observed at all densities and changes 
systematically with the peak electron density. In Fig. 3 
we show several examples. The decrease in the slope with 
time is more pronounced at higher densities. Meanwhile, 
the simultaneously obtained expansions of the electron 
density profile (not shown) are all linear and similar to 
Fig. 2c. To approximately compare the rates of spin dif- 
fusion and carrier diffusion, we perform linear fits to the 
data with t > 150 ps (solid lines in Fig. 3). In this range 
the expansion is approximately linear. The spin diffusion 
coefficients deduced by the fits are plotted in Fig. 4 (cir- 
cles) as a function of peak electron density. The ambipo- 
lar carrier diffusion coefficients are also plotted (squares) 
for comparison. In the density range studied, the am- 
bipolar carrier diffusion coefficient is almost constant. In 
contrast, the spin diffusion coefficient decreases signifi- 
cantly with peak electron density. 
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FIG. 3: (Color online) Expansions of spin density profiles at 
several peak electron densities as indicated in each panel, ob- 
tained by repeating the procedure summarized in Fig. 2. The 
sublinear expansion is generally observed. Linear fits to data 
with t > 150 ps (solid lines) are used in order to approxi- 
mately illustrate the density dependence of spin diffusion. 



IV. DRIFT-DIFFUSION THEORY 

We denote by p the density of holes. The drift-diffusion 
equations for the three coupled densities, n.f, uj,, and p, 
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where Cj, a± and a p are the ordinary conductivities of 
spin-up electrons, spin-down electrons and holes respec- 
tively; D aa i is the spin-diffusion matrix for electrons; E 
is the electric field; l/r-j-; and 1/rj.f are the spin flip rates 
from up to down and from down to up, respectively; r r 
is the electron- hole recombination time. 

Some physical assumptions underlie the above equa- 
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FIG. 4: (Color online) Ambipolar carrier diffusion coefficient 
(squares) and spin diffusion coefficient (circles) as functions 
of peak electron density obtained by repeating the procedure 
summarized in Fig. 2 at various electron densities. These 
coefficients are deduced by linear fits as shown in Figs. 2 and 
3. Significant decrease in spin diffusion coefficient with peak 
electron density is observed. The carrier diffusion coefficient 
is almost constant. 



tions. First of all, we have completely neglected the spin 
polarization of the holes, assuming that the spin relax- 
ation time of the holes is very short. Second, we have 
ignored the influence of carrier thermalization and relax- 
ation on the diffusion process, and assumed a constant 
carrier temperature. This is justified since these pro- 
cesses are expected to end on the time scale of 10 ps even 
with the influence of hot-phonon effect.— Indeed, the lin- 
ear expansion of carrier density profile observed in Fig. 2c 
confirms that the diffusion coefficient is a constant, and 
the influence of thermalization and relaxation on diffu- 
sion is negligible although the electrons are initially ex- 
cited with an excess energy of 40 meV. Finally, by keep- 
ing the off-diagonal elements of the spin-diffusion matrix 
(e.g. D^i) we have allowed in principle for the effect of 
spin Coulomb drag,— whereby a gradient of spin-up den- 
sity can drive a current of spin down and vice- versa. The 
explicit form of this matrix is 



D = 



1 + 7T 



1 i '*T "T 

1 H L 7T —t-^VT 

n i n ' 

n l 1 i i. 

— 7T 1 H L 7T 



(5) 



where D n = finksT/e is the electron diffusion constant, 
t is the momentum relaxation time for electrons, and 7 is 
the spin Coulomb drag coefficient. The steps for deriving 
the spin diffusion matrix from the spin resistivity matrix 
are outlined in the Appendix. 

The equations for n-t and n± can be rewritten in terms 
of total density n = n-\ + 714 and spin density s = n-j- — 7i| 



as follows: 

dn n 

— =V • [/i n nE + D n Vn] 

at tv 



(6) 



^ = V • [fi n sE + D s ( Vs + sjtV Inn)]- — -—, 
at t s r r 

(7) 

where /Lt„ and \x v are the mobilities of electrons and holes, 
respectively, and we have used the standard relations 
a t(l) = n t(l) e t 1 n and cr p — epfip. In the second equa- 
tion we have set 



D. = 



Dn 

1 + 7T 



(8) 



and have assumed = I/t^j. Finally, + = 

1/t s , where l/r s is the spin relaxation rate. 

The electric field E arises from the small charge im- 
balance that inevitably occurs as low-mobility holes try 
to keep up with high-mobility electrons. Even though 
the charge imbalance is very small, it would not be legit- 
imate to neglect this field, since the conductivity of the 
electrons is very high. Comparing the equations for n 
[B] and p [3] , and making use of the approximate charge 
neutrality condition p ~ n, we identify the electric field 
as 



E 



D p - D n Vn 



(9) 



Given the relative magnitudes of D p < D ni we can see 
that the electric field points away from high density re- 
gions. The impact of this small electric field is a slight 
resistance to the diffusion of electrons. 

Substituting this into Eqs. ([6]) and ([7]) we arrive at our 
final equations: 



dn _ ~ n 
— =D a V 2 n - - 
at tv 

% =V- [(A,- D 8 )sV Inn- 
at 



n Vsl - - - - 



where 



D a 



l^n H~ ftp 



(10) 

(11) 

(12) 



is the ambipolar diffusion constant, which is intermedi- 
ate between D n and D p but numerically closer to the 
diffusion constant of the less mobile species (holes in 
this case). Equation (fT0|) is the standard electron diffu- 
sion equation with ambipolar diffusion constant. Strictly 
speaking, the recombination time can vary with density, 
complicating the solution. Experimentally, a recombina- 
tion time between 160 ps at low densities (~ 10 17 cm~ 3 ) 
and 190 ps at higher densities (~ 6 x 10 17 cm~ 3 ) was 
observed. For the time being, consider the recombina- 
tion time to be spatially uniform, a point to which we 
will return and justify later. Assuming that the electron 
density packet at the initial time t = has a Gaussian 
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shape, n(r, 0) = iVe 41n ( 2 ) r / w o : we see that the solution 
to Eq. (fTU|) at time t is given by 

„( r>t ) = JV ("^Ly /2 e-41n(2)^/^( t )- t /r r; (lg) 

where w^(t) = Wq + 161n(2)D a i. The linear growth in 
time of the area covered by the packet is in excellent 
agreement with the experimental observation, as shown 
in Fig. 2c. 

Let us now consider Eq. pip for the spin density. We 
assume that the initial spin distribution is proportional 
to the density distribution, i.e., we have 

s(r,0) = C*n(r,0) (14) 

where C is a constant independent of position. For the 
present experiment the expected value of C is C = 1/2. 
The number of spin-up electrons is three times larger 
than the number of spin-down electrons, owing to selec- 
tion rules for two-photon absorption of circularly polar- 
ized light. 42 Then, if the spin relaxation time is neglected 
we see immediately - by direct substitution - that the 
solution of Eq. (JTTJ) is 

s{r,t) = Cn(r,t) . (15) 

In other words, the spin density and the ordinary density 
diffuse at exactly the same rate, controlled by D a . 

It should be noted that this result holds irrespective of 
the value of the spin Coulomb drag coefficient. The am- 
bipolar spin diffusion in our configuration is in sharp con- 
trast with the case of unipolar spin packets in spin-grating 
experimentsi 25 ' 26 i 27 ' 28 i 29 In those experiments, two laser 
pulses that are oppositely linearly polarized are used to 
excited carriers. Since the interference does not intro- 
duce any intensity modulation of the light, the total elec- 
tron density and hole density excited are not modulated. 
Hence, there is no density gradient to cause diffusion of 
holes and electrons. However, the interference gives rise 
to spatial modulations of the light polarization, resulting 
in periodic variations in spin density of electrons. The ex- 
cited spin grating decays due to the diffusion of spin-up 
and spin-down electrons from high density to low density 
regions, within the laser spot. In this case, the same the- 
ory predicts the spin diffusion constant to be D s , given 
by Eq. ([8])i=a In this unipolar spin diffusion regime, the 
requirement of charge neutrality forces the electrons of 
opposite spin orientations to move in opposite directions: 
spin Coulomb drag arises from this relative motion. In 
the ambipolar spin- diffusion process, the spin-density 
profile overlaps with the hole density profile, and charge 
neutrality is ensured by the holes. Thus, the electrons 
move in the same direction regardless of spin - there is 
no relative motion and therefore no spin Coulomb drag, 
and the diffusion is controlled by the ambipolar diffusion 
coefficient. 

The above solution (|15p was found for the initial con- 
dition (|14p but we have checked that, even if the spin 



density is not initially proportional to the density, it even- 
tually becomes proportional to the density at sufficiently 
long times. 

We conclude that the experimental observation of an 
apparently decreasing spin diffusion coefficient cannot be 
explained within the framework of the drift-diffusion the- 
ory, unless one is willing to include the spin relaxation 
time. Further, it is easy to see that a homogeneous spin 
relaxation time (l/r s independent of density and hence 
of position) will not change the situation, for the solution 
of Eq. (JTTJ) in the presence of such a relaxation time is 
simply 

s{r,t) = Cn(r,t)e- t/T ° . (16) 

There is now decay of the spin due to spin relaxation in 
addition to the global decay of carrier and spin densities 
due to electron-hole recombination. 



V. INHOMOGENEOUS SPIN RELAXATION 

In view of the above discussion, we now examine the 
possibility of explaining the experimental data in terms of 
non-homogeneous spin relaxation. Suppose, for instance, 
that the spin relaxation rate were larger at lower density, 
i.e., larger in the tails of the spin packet than at its cen- 
ter. This would produce the impression of slower spin 
diffusion, since the outward diffusion of the spin would 
be hidden by the more rapid decay of the spin at the 
edges (Fig. 5). To verify this idea, we measure the spin 
relaxation time as a function of peak electron density 
by using the same pump-probe technique. The results 
shown in Fig. 7 confirm this conjecture (see discussion 
below). The spin relaxation time is larger for packets of 
higher density, suggesting that l/r s increases with de- 
creasing density. However, these are measurements of 
the lifetime of the spin packet as a whole, whereas in the 
drift-diffusion theory we need a position-dependent spin 
relaxation rate, determined by the local density. 

We have developed a model for the position depen- 
dence of the spin relaxation rate along the following 
lines. First, we notice that electron- impurity effects 
or spin-orbit interactions with the lattice could not ac- 
count for the spatial variation of l/r s within the packet, 
since the impurity environment and the crystalline en- 
vironment are uniform over the region occupied by the 
packet. What is not uniform is the frequency of electron- 
electron (e-e) and electron-hole (e-h) collisions. Because 
the Coulomb interaction is screened and thus effectively 
becomes of short range, we can say that the carriers near 
the center of the packet experience an electronic envi- 
ronment of higher density, and their scattering rate is 
accordingly higher than in the tails of the packet, where 
they experience an electronic environment of lower den- 
sity. 

How does the variation in the carrier-carrier scattering 
time translate into the observed spatial variation in the 
spin relaxation rate? A plausible answer comes from the 
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Distance from center of packet 

FIG. 5: (Color online) As electron density decreases, the spin 
relaxation rate increases (dashed line) . A non-uniform spin re- 
laxation rate across the packet gives the appearance of slower 
width increase. The length of the arrows is given by the prod- 
uct of density and relaxation rate, (linear l/r s for illustrative 
purposes only) 



D'yakonov-Perel' spin relaxation mechanisi m 12 ' In this 
mechanism, the spin relaxation is primarily due to the 
spin-orbit interaction with the lattice. In GaAs the pri- 
mary mechanism is the independent precession of each 
electron in the Dresselhaus fc-dependent effective mag- 
netic field^ 

n k =a c K 2 (2m* 3 E g ^ k, (17) 
k = {k x (ky — fe z ) , ky [k 2 — k 2 ^ , k z (fc 2 , — fc^)} (18) 

where the parameter a c ~ 0.07 for GaAs, m* is the effec- 
tive mass of electrons in the conduction band, and E g is 
the band-gap energy. Notice that the effective magnetic 
field depends on the momentum k of the electron, e-e and 
e-h collisions change the momentum of the electron, and 
thus change the direction and magnitude of fl^. This 
hinders the process of spin relaxation, especially in re- 
gions where momentum-changing collisions are frequent, 
i.e. high-density regions. A standard analysis for cubic 
systems leads to the conclusion that 

1 39 e 3 

- = <n Xfe r )^— Tac — (19) 

where r* is related to the momentum relaxation time 
by the effectiveness of the scattering event on the ran- 
domization of the axis of fifc, £l±k is the component of 
Ctk perpendicular to the spin axis, ek — h 2 k 2 /2m* is the 
one-particle energy, and the angular bracket denotes a 
thermal average in momentum space. 

We have evaluated the contribution of e-e and e-h in- 
teractions to the lifetime of a momentum state in a homo- 
geneous environment of density n according to standard 
formulas reported for example in Ref. 56 , but evaluated at 
the parameters of the experiment. In order to simplify 
evaluation of the momentum relaxation time we use a 
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FIG. 6: (Color online) The momentum relaxation times due 
to e-e (solid line) and e-h (dashed line) scatterings are on the 
same order of magnitude as the background scattering (dotted 
line) in high density regions (> 10 17 cm -3 ) but become quite 
long at lower densities. Under the present interaction model, 
the rate of carrier-carrier collisions is proportional to n. Here 
the relaxation time is calculated for wave vector k — 1.43 x 10 s 
m -1 and screening wave vector ko = 1.89 x 10 8 m _1 . 

screened Coulomb interaction model of the form 

p 2 

— e - kor 
e r r 

where ko is a fitting parameter of the order of the De- 
bye screening wave vector at n = 10 17 cm -3 and Boltz- 
mann statistics are used throughout. This is admittedly 
a significant approximation as the electrons are degener- 
ate above ~ 10 17 cm -3 . While numerical accuracy may 
be improved, qualitative aspects are generalizable. Our 
results are shown in Fig. 6, where the momentum relax- 
ation times due to e-e and e-h scatterings are compared 
to the background scattering due to impurities, which is 
density independent. These values have been obtained 
for a typical thermal value of k — 1.43 x 10 s m -1 with a 
screening wave vector of ko = 1.89 x 10 8 m _1 . Applying 
Matthicsscn's rule, high density regions exhibit a higher 
rate of collisions which could potentially slow spin relax- 
ation, whereas low density regions are governed primarily 
by lattice scattering. 

Using Eq. (flT))) , we arrive at t s shown in Fig. 7, where 
the calculated spin relaxation time (solid line) is shown 
to be in good agreement with experimentally obtained 
t s (circles). We remind the reader that experimental t s 
is plotted as a function of peak electron density, while 
Eq. (fTTJ|) is being applied on a microscopic scale. At this 
point, the earlier approximation of a spatially uniform 
tv can be justified. Whereas the range in spin relax- 
ation time for 10 17 cm~ 3 < n < 6 x 10 17 cm -3 varies by 
as much as 600 ps, the electron-hole recombination time 
only varies by 30 ps. Clearly, the spatial dependence 
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FIG. 7: (Color online) The spin relaxation time of electrons is 
inversely proportional to the momentum relaxation time, thus 
increases with density. Experimentally measured t s (circles) 
are plotted versus the peak electron density. 



of the recombination rate is swamped off by the spatial 
dependence of the spin relaxation. 

The results of the calculation give us confidence that 
the D'yakonov-Perel' mechanism, limited by the carrier- 
carrier scattering, can indeed produce a spin relaxation of 
the right order of magnitude and - more importantly for 
us - one that is faster in the tails of the packet. We note 
that our results are consistent with a recent study where 
electron spin relaxation time in bulk GaAs at room tem- 
perature was found to increase linearly with carrier den- 
sity when the D'yakonov-Perel' mechanism dominates. 57 

With the inhomogeneous spin relaxation time, the spa- 
tiotemporal evolution of a spin polarized packet was eval- 
uated using Eq. with parameters similar to those 
used in Fig. 2. One dimensional dynamics were used for 
the simple goal of comparing diffusion rates. Figure 8 
illustrates the squared widths of the profiles of electron 
density and spin density over an interval of 300 ps, calcu- 
lated by fitting a Gaussian curve at several times during 
the evolution. An interesting consequence of the inho- 
mogeneous spin relaxation is that by breaking the triv- 
ial solution it "turns on" the spin Coulomb drag, which 
would otherwise be completely inoperative. The spin dif- 
fusion rates with and without spin Coulomb drag are also 
compared. The differences are rather small, reflecting the 
fact that the slowing down of the spin is mostly apparent, 
i.e. due to the loss of spin in the edges rather than to the 
relative motion between up and down spin components. 
Nonetheless, a small effect is visible and goes in the di- 
rection of further reducing the apparent spin diffusion. 




50 100 150 200 250 300 



Time, t (ps) 

FIG. 8: (Color online) Squared width of the spin density 
packet calculated with (dotted line) and without (dashed 
line) taking into consideration the spin Coulomb drag. These 
are considerably smaller than the squared width of a density 
packet (solid line). It is evident that the spin Coulomb drag 
plays only a minor role in further reducing the apparent diffu- 
sion rate. Parameters used matched those of the experiment 
in Fig. 2. 



VI. SUMMARY 

We have studied, both theoretically and experimen- 
tally, the diffusion of optically injected spin-polarized car- 
riers in undoped GaAs quantum wells at 80 K. The exper- 
iment is performed with a high-resolution optical pump- 
probe technique. Spatiotemporal dynamics of locally in- 
jected spin-polarized carriers are directly resolved. By 
comparing expansions of electron and spin density pro- 
files, we found that the spin diffusion cannot be described 
as a classical diffusion process with a constant diffusion 
coefficient. The spin diffusion appears to be slower than 
the ambipolar carrier diffusion. Our theoretical analy- 
sis, based on a three-component drift-diffusion equation, 
shows that the long-time behavior of the spin density 
can be understood in terms of a spin-relaxation rate that 
grows with decreasing density. This behavior is consis- 
tent with a model of D'yakonov-Perel' relaxation limited 
by the Coulomb scattering between carriers. However, 
the short-time behavior of the time varying spin diffu- 
sion coefficient (within 100 ps) remains for the time be- 
ing beyond the reach of our drift-diffusion theory, and is 
subject to our further investigation. 
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APPENDIX A: SPIN DIFFUSION MATRIX 

The spin diffusion matrix is derived from the resistivity 
matrix by the application of Einstein's relations. The 
resistivity matrix is given by^ 



— H — -T7 

— T7 



— T7 

n ' 



(Al) 



where we have assumed no external electric field. In ad- 
dition, the momentum relaxation rate due to electron- 
impurity collisions which flip the spin was taken as negli- 
gible in comparison to non-flip processes, r is essentially 
the Drude scattering time. 

Inverting the resistivity to find conductivity (tr), Ein- 
stein's relation gives the diffusion matrix 



e 2 D a p = yv 

7 



[x 



I 7/3 



(A2) 



where the spin susceptibility matrix \ is approximated 
as 



(A3) 
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The diffusion matrix is finally. 
D 



D = 



1 , nt n-r 
1 + n T T n T "f 

— T7 1 H L T r y 

n 1 n 



1+JT 

where the following relation was used, 



Mr, 



(A4) 
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